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NOMENCLATURE 
a Ideality factor of the diode 
G Solar irradiation 
Gy Solar irradiation at STC condition 
f1 Output current of PV cell 
la Shockley diode current 
Linp Current of the cell at MPP 
i Reverse saturation current of the cell 
lon Nominal reverse saturation current of the cell at STC condition 
lov Incident light current of the cell 
lovn Nominal incident light current of the cell at STC condition 
ee Short circuit current of the cell 
lees Nominal short circuit current of the cell at STC condition 
k Boltzmann constant (k=1.3806503*10723 J/K) 
KCL Kirchhoff’s current law (KCL) 
K; Temperature coefficient of I,, 
K, Temperature coefficient of V,. 
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Electron charge (1.60217646*107?? C) 
Series resistance 

Previous value of R, 

Parallel resistance 

Previous value of R,, 

Initial value of R,, 


p,min 
STC Standard Test Condition 
T Absolute temperature of the cell in Kelvin 
Te Ambient temperature. 
L Absolute temperature of the cell in Kelvin at STC condition 
le| Absolute error 
V Output voltage of PV 
Vinp Voltage at MPP 
Vee Open circuit voltage 
Vea Open circuit voltage 
AP es Difference between Prax e¢ and Praxm 
APmaxpre The previous value of APygx 
Vi Thermal voltage at 25 C° 
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Np Number of parallel cells 
Ns Number of series cells 
ieee Maximum power point provided by datasheet 
Pe ae Maximum power point provided by model 


1, INTRODUCTION 

A remarkable percentage of the used energy is produced by the fossil fuel. The COz and many other 
harmful gases that produced as a result of the fuel's combustion cause many problems to the earth's environment 
such as Global Warming. Meanwhile, some studies prove that the fossil fuel has an end, therefore, the earth 
needs alternative energy resources. Sun is one of the most remarkable and renewable energy resources which 
can produce a useable clean energy. The Sun gives us (approximately) 170 PW [1]. Solar energy might not be 
the best solution nowadays because of the high cost but it might be for our future generation. Furthermore, 
a study shows that the renewable resources are growing fast and they become the second largest sources in 
2015 and the first energy resources in 2035 [2]. Solar energy is the fastest growing source among them and it 
doubles each two years [3]. 

Solar cells convert the light from the sun into DC electrical energy and this occurs in some materials 
which have the ability to capture photons and release electrons and then current flows. Silicon is the main 
material in the industry to produce photovoltaic (PV) modules. In this paper, mathematical modelling 1s used 
in order to analyse the performance characteristics of the PV modules. Using the mathematical modelling 
makes the analysis much easier in contrast to the physical modelling [4]. By making mathematical modelling 
for the PV modules, manipulating and changing the parameters of a given module are easily achieved [5]. 

The output values (voltage, current and power) of PV models depend on irradiation, temperature 
and load current. The variations of these (irradiation, temperature and load current) play an important role in 
the modelling [6]. The single-diode model is not the only available model. In some papers, two and three-diode 
models are presented. The two-diode model in [7-11] 1s used for having higher accuracy. The second diode is 
added to involve the effect of the carrier recombination. The third diode is added in the three-diode model in 
[12] to involve the missing effects in the two-diode model. Single-diode model [13] is the most popular model 
because of the simplicity. It is categorized into two types which are four-parameter and five-parameter models. 
These five parameters including the parallel resistance are series resistance, incident light current, diode reverse 
saturation current, and ideality factor. The five-parameter model is considered more accurate and widely 
applied than the four-parameter model [14]. 

To extract these five parameters, several methods were applied. One method is to identify 
the parameters using the curve-fitting for the experimental data collected under specially controlled 
environment (SCE) such as sun simulator [8]. Deferent methods such as neural network (NN) 
and neuro-fuzzy model are used. NN model [15, 16] has a drawback which is the need of much data 
for training the neural network while neuro-fuzzy model needs fewer data [17]. 

Another method is to extract the parameters from driven equations with the help of the electrical 
specification data of the solar module provided in datasheets [18, 19], which are: the nominal short-circuit 
current (U,.,), nominal open-circuit voltage (V.,), experimental peak output power (Pygx.), voltage 
at MPP (Viny»), current at MPP (/,,,,), the temperature coefficient of the open-circuit voltage (K,), 
the temperature coefficient of the short-circuit current (K;) and the temperature at the nominal operating cell 
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temperature (NOCT) under Standard Test Condition (STC). This method is much more reliable because 
the parameters of the model are computed mathematically using reliable data (datasheet). 
Therefore, it produces more accurate models, but in the price of time. It uses iteration methods, therefore, 
it consume too much time to find the model parameters. This paper tries to solve the unbearable 
time-consumption. In the meantime, it produces an accurate model. In this paper, we start with 
the methodology that includes model of a solar cell and the proposed modelling algorithm (section 2). 
Then, KC200GT solar array is used in section 3 to examine the EFPPM’s accuracy and speed and compare 
with Villalva’s modelling algorithm [24]. Lastly, section 4 displays the conclusion. 


2. RESEARCH METHOD 
2.1. Model of a solar cell 

There are two popular types of cell models. The first model depends on one diode and it is called 
single-diode model [20, 21]. While the second type depends on two diodes and it is called double-diode 
model [8]. However, the double-diode model has better accurate results, but it has too many complicated 
equations. The single-diode model has also accurate results and easier mathematical expressions [22]. 
Therefore, Villalva’s algorithm is based on single diode model as in Figure 1. 





Figure 1. Practical solar cell model 


Villalva’s algorithm has many improvements for calculating the PV model parameters as will be 
clarified later. Therefore an accurate and reliable PV model can be obtained. This method depends mainly on 
finding the model parameters using three point in the I-V curve. The three points which are open circuit, short 
circuit and maximum power point are easily obtained by the data provided in datasheets under STCs. The I-V 
equation for a PV array 1s: 


I = Ipy — I, [exp (= *) - 1] ae (1) 


where [,,, and I, are the photocurrent and the diode reverse saturation current, respectively, q is the electron 
charge (1.6 x 1071? C), k is the Boltzmann constant (1.38 x 10723 J/K), a is the modified ideal factor, R, 
and R, are, respectively, the series and parallel resistance (Q), and T is the cell temperature (K). The (1) 
is obtained from Figure 1 based on Kirchhoff's current law (KCL), where the diode current [, 
is substituted by; 





- QqV+IRs _ 
la = lo lex (Se) — 1 @) 
Adjusting (1) at the three points can lead to find the unknown parameters [,,,, [,, Rp, Rs, and a. Ip, can be 
found by (3). 
G 
Lg — Upvn =K;7T —T,)| Gy (3) 


It is generally assumed [,,=/,, (Short circuit current of the cell) because the R, is low and R,, is high. 
The diode saturation current 1s normally expressed as [23]. 
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lo= lon () em [EP (=-- 7) 4) 


The (4) can be replaced by (5) and provides an improvement for the PV model because the dependency 
on temperature is included using the voltage and current temperature coefficient (K,, and K;) [24]. 


(Voen+KyAT ) / | 
av; —1 


lL, = 
exp| 


(5) 


Finding R,, and R, is done based on that there is a point where the model maximum power point Fraxm equals 
the datasheet maximum power point Pygye¢ as in (6) [24]. Adjusting the I-V equation by incrementing R, starts 
from zero and then finds R,, by (7). Incrementing R, should be linearly with a small and fixed step increase in 
order to get the desired values of R, and R,. 


(VinvtlmpRs) @ VinptlmpRs 
P — V f = | Jex [ee s) a) — | — Yolo — Pp (6) 
max,m mp pV O p aN, kT Ry max,e 
(Vinp+Imp Rs) 
R= Pr (7) 
aay (VmptimpRs) q 4 ,—Pmax,e 
pv QO p aNs kT o Vmp 








Further improve is developed by [24] using (8) which warrants more accurate matching between model I-V 
curve and real PV array I-V curve so I, # Isc. 


Rp +Rs 
Lown > Rp Len (8) 





The initial value of R, is O while the intial value of R,, is obtained by: 





Vinp Vocon—Vmp 
Kpimin 7 Iscn—lmp - “ty ”) 
2.2. Proposed algorithm 

EFPPM is designed similar to Villalva’s algorithm. However, in the designed EFPPM, dynamic 

incrementation is proposed instead of the static incrementation proposed in Villalva’s algorithm. 
The dynamic incrementation is used to improve the speed of algorithm and reduces the computation time when 
finding the model parameters. The initial incrementation step in EFPPM is set to 0.1 (0.1 1s considered a big 
step value because the value of R, 1s small). Incrementing continues until the value of R, jumps over the desired 
value. Then the last value of R,; before exceeding the desired value will be increased using a new increment 
step (the new increment step=the previous increment step divided by 10). These processes continue until 
matching the Pygym With Prax. The ideas are to use a dynamic increment step for R, incrementation starting 
from a big value 0.1 until the best value of R, and R,, are found. The static increment step is replaced by the 
dynamic increment step to reduce the iteration numbers of the modelling algorithm in [24]. Figure 2 shows 
the flowchart of the proposed EFPPM. 


3. RESULTS AND DISCUSSION 

EFPPM reduces the number of the iteration and has the same accuracy of the Villalva’s algorithm. 
Reduction is done based on using dynamic incrementation. Table 1 shows the electrical characteristics 
of KC200GT solar array [25] provided by the datasheet. The KC200GT solar array is used here to create 
a PV model using EFPPM and the Villalva’s algorithm and then makes a comparison between them. 

Figures 3 and 4 show the I-V and P-V curves plotted for different values of R, and R, using both 
the Villalva’s algorithm and EFPPM, respectively. As we notice, the plotted curves in Figures 3. (a) and (b) 
move toward left as R, is slowly increased. In Villalva’s algorithm, R, starts from zero. The desired R, and R,, 
are found when Pygym matches with Pyg,. using equation (6) as shown in Figure 3. (b). In Figures 4. (a) 
and (b), the plotted curves swing around MPP until the Pyg,.m matches with Pyg,.. Table 2 shows how R, is 
changed in EFPPM. The desired R, and R, are 0.221 Q and 415.76 Q, respectively. Both methods 
the Villalva’s algorithm and EFPPM produce same result of the desired R, and R, as shown in Figures 5. (a) 
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and (b) but with different significant speed. They show that the model curves I-V and P-V exactly match with 
the experimental curves I-V and P-V of KC200GT solar array given by the datasheet at three points: short 
circuit point, MPP and open circuit point. 





Inputs:T, G,I,,eq(5) 
R, =0,APnay = lel +1 
Rs pre» Ry ore = 0 
Increment = 0.1 
Ry = Ry min eq(9) 
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Increment R,; = R, + Incremnt 











Figure 2. Flowchart of EFPPM 


Table 1. Typical electrical characteristics of KC200GT module under STCs 


parameter value 
Paves 200.143 W 
Vinp 26.3 V 
Lag 7.61 A 
ee 8.21 A 
Voen 32.9 V 
K; 0.0032 A/K 
K, -0.1230 V/K 
Ns 54 
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Figure 3. (a) I-V and (b) P-V curves plotted for many values of R, and R,, for the Villalva’s algorithm 
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Figure 4. (a) I-V and (b) P-V curves plotted for many values of R, and R, for EFPPM 
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Figure 5. (a) I-V and (b) P-V curves plotted for the optimum value of R, and R,, for both EFPPM and the 
Villalva’s algorithm. 


The number of plotted I-V and P-V curves until finding the optimum values of R, and R, using 
EFPPM is 10 while using the Villalva’s algorithm is 221. In terms of computation time, EFPPM requires 5.7 
Sec. when executed in the processor (2.0 GHz Intel Core 2 Duo T5800 compared to 123.4 sec. for Villalva’s 
algorithm. This much time comes due to solving the nonlinear equation (1) for I € [0, I all and V € [0, Vonk 
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I5en and Vion are 8.21 and 32.9, respectively. The values of R, used in the iteration process using EFPPM 
are listed in Table 2. 


Table 2. R, Incrementation values of EFPPM 


R, Values 
First Sequence Second Sequence Third Sequence 
0 0.20 0.220 
0.1 0.21 0.221 
0.2 0.22 


The desired parameters values with a comparison between EFPPM and the Villalva’s algorithm in 
term of iteration no and computation time are listed in Table 3. The table shows that both methods give same 
result but different computation time. Table 3 shows the superiority for EFPPM against Villalva’s algorithm 
in term of time efficiency. Lastly, despite of the high accuracy of Villalva’s algorithm, it consumes too much 
time to find the five parameters of the PV modules. Therefore, we need a modelling algorithm that produces 
high accurate PV models with acceptable computation time and that is available in EFPPM. EFPPM can be 
more effective with large solar arrays that have large R,, because the reduction of time computation done 
by EFPPM will be more noticeable. 


Table 3. Parameters of the adjusted model of the Kc200gt solar array at nominal operating conditions 


Parameters Villalva’s algorithm [24] EFPPM 
Pied 200.143 W 200.143 W 
Vinp 26.3 V 26.3 V 
Inp 7.61 A 7.61 A 
fea 8.21 A 8.21 A 
Veen 32.9 V 32.9 V 
a 1.3 iW) 
sie 8.214A 8.214 A 
lay 9.8225e-08 9.8225e-08 
R, 0.221 0.221 
Ry 415.78 415.78 
Iteration Num 221 10 
Computation Time 123.447 Seconds 5.725 Seconds 


4. CONCLUSION 

Solar energy has some drawbacks such as high cost, recombination, shading, and day sunlight. Many 
researchers try to maximize the efficiency of photovoltaic system by studying the I-V and P-V characteristics 
curves. Mathematical modelling gives us an easy and accurate tool to analyse I-V and P-V characteristics 
curves. In this paper, the main objective is to overcome the slowness of the existing modelling algorithms. The 
proposed modelling algorithm (EFPPM) is a method to find the five unknown PV model parameters (a, 
I,, ny, Rs, and R,,). Finding the PV model parameters is done based on adjusting the I-V curve in three points 
that are provided by datasheet (short circuit, open circuit, and MPP). An iterative process of increasing R, and 
then find R,, and I, 1s the main key to adjust the I-V curve. This iteration processes, based on fixed increment, 
consumes much time because of solving nonlinear equation. An improvement is fulfilled by minimizing 
the number of the iterations. Instead of using fixed increment to find the optimum values of PV model 
parameters, a dynamic increment is used. The time consumed to find the model parameters of KC200GT solar 
array 1s 5.725 seconds using EFPPM instead of 123.447 seconds using a similar work. This result shows that 
EFPPM has superiority in term of the computation time. 
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